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SUMMARY 


Accurate  prediction  of  the  lifetimes  of  satellites  decaying  naturally  under 
the  action  of  air  drao  is  one  of  the  most  difficult  problems  of  orbital  dynamics. 
Lifetime  estimation  is  also  a  matter  of  practical  importance  when  a  satellite 
perceived  as  'dangerous'  is  about  to  decay.  The  time  when  the  satellite  finally 
plunges  into  the  lower  atmosphere  is  affected  by  at  least-height*  identifiable  types 
of  variation  .n  upper-atmosphere  density,  by  odd  zonal  harmonics  in  the  gravita¬ 
tional  field,,  and  by  possible  changes  in  the  mode  of  rotation  of  the  satellite. 

For  use  in  the  RAE  Table  of  Satellites,  we  have -since  1958, made  (and  revised) 
lifetime  estimates  for  all  satellites,  rockets  and  other  major  components  (but 
not  fragments),  More  tnar.  10000  such  predictions  have  been  made,  nearly  all  being 
calculated  from  the  current  observed  decay  rate.  This  Report  summarizes  the 
methods  used  in  the  past  and  gives  the  results  of  new  orbital  theory  developed  to 
take  better  account  of  the  effects  of  atmospheric  oblatencss  and  odd  zonal 
harmonics  in  the  aeopotential.-^  Hie  Report  also  describes  new  approximate  methods 
to  allow  for  the  effect  of  the  solar-cycle  and  day-to-night  variations  in  air 
density,  and  includes  discussion  of  other  sources  of  error.  There  are  separate 
sections ’giving^ techniques .for  use  when  no  observed  values  of  decay  rate  are 
available.  'Numerical-integration  methods  are  not  covered, 'but* ye  give  the 
graphical-analytical  methods  devised  for  two  common  types  of  high-eccentricity 
orbit  -  geosynchronous  transfer  orbits  and  Molniya-type  orbits.  — 
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1  INTRODUCTION 

Since  1958  the  likely  lifetimes  of  all  spacecraft  and  rockets  in  orbit  have 

been  estimated  and  updated  at  RAE  for  use  in  the  RAE  Table  of  Earth  Satellites* . 

.  2 

For  orbits  of  low  eccentricity,  the  methods  have  been  based  on  the  current 
observed  decay  rate  and  orbital  theory^  derived  on  the  assumption  that  the  air 
density  has  no  time  variations;  but  a  correction  has  been  made  for  the  semi" 

4 

annual  variation  m  air  density.  Recently  ,  the  theory  has  been  extended  to 
derive  corrections  for  the  oscillations  m  perigee  height  (of  order  10  km)  caused 
by  the  odd  zonal  harmonics  m  the  geopotential  and  by  atmospheric  oblateness. 

The  present  Report  incorporates  the  corrections  for  these  two  effects  and  also 
derives  approximate  correction  factors  for  the  effects  of  long-term  solar  activity 
and  the  day-to-night  variation  in  density,  giving  a  total  of  five  corrections  to 
the  basic  lifetime  estimate.  The  effects  of  short-term  variations  in  air  density 
are  still  unmodelled,  however,  because  they  are  unpredictable,  and  they  may  be 
expected  to  produce  errors  of  up  to  I0Z  in  the  prediction  of  the  remaining  life¬ 
time  of  the  satellite,  unless  solar  activity  remains  steady. 

For  high-eccentricity  orbits,  lunisolar  perturbations  need  to  be  taken  into 
account.  Although  numerical  integration  is  usually  needed  for  these  orbits,  two 
analytical  methods  have  been  developed  to  estimate  lifetimes  approximately,  for 
low-inclination  transfer  orbits'*  and  for  Molniya-type  orbits**.  For  completeness, 
.he  necessary  graphs  for  using  these  methods  have  also  beer  included  here,,  with 
some  small  improvement  on  the  original  formulations. 

Mere  accurate  predictions  for  high-eccentric^ty  orb. is  can  be  made  by 
numerical  integration  of  the  equations  of  motion,  although  this  integration  often 
has  to  continue  for  more  than  25  years  and  is  therefore  very  demanding  in  com¬ 
puter  time.  For  these  satellites  the  zero-drag  PROD  computer  program^  is  used 
for  making  the  first  lifetime  estimate,  which  is  often  quite  accurate  when  perigee 
rises  to  a  great  height  and  then  cleanly  descends  to  below  100  km,  but  is  less 
satisfactory  if  the  perigee  tends  to  oscillate  at  heights  of  150-250  km  As 
decay  approaches,  we  use  the  program  PTDEC®,  which  includes  drag  via  an  atmos¬ 
pheric  model  but  is  much  more  expensive  in  computer  time.  The  use  of  these 
9 

programs  is  described  elsewhere  ,  and  they  are  not  discussed  here. 

It  is  of  course  possible  to  predict  the  lifetime  of  a  low-eccentricity 
orbit  by  numerical  integration  instead  of  the  graphical  methods  presented  here. 

But  the  program  for  numerical  integration  must  incorporate  an  atmospheric  model 
that  covers  as  many  variations  as  the  graphical  method  if  tne  results  are  to  be 
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as  accurate  s  those  of  the  graphical  method.  (Unfortunately  PTDEC  does  not 

include  semi-annual  or  day-to-night  effects,  or  variations  in  solar  activity, 

9 

and,  even  so,  cannot  m  practice  be  run  for  long  enough  .) 

It  would  also  be  possible  to  computerize  the  graphs  given  here:  but  the 
calculations  for  lifetime  estimation  are  trivial  with  the  graphical  method,  once 
the  decay  rate  is  known.  The  'difficult'  parts  of  the  task  are  deriving  a 
reliable  decay  rate,  and  arranging  to  be  able  to  lay  hands  quickly  on  several 
sets  of  orbital  elements  for  any  of  the  6000  satellites  in  orbit.  The  methods 
presented  here  allow  lifetimes  to  be  predicted  without  computer  capability. 

The  arrangement  of  the  Report  is  as  follows:  after  a  digression  m 
section  2  on  the  'threat'  of  satellite  decay,  section  3  gives  the  basic  method 
for  low-eccentricity  orbits,  in  terms  of  the  currently-observed  decay  rate;  then 
sections  4  to  7  detail  the  corrections  required,  which  are  summarized  in  section  8, 
where  other  sources  of  error  are  also  discussed.  Section  9  gives  the  alternative 
method  to  be  used  when  values  of  the  current  decay  rate  are  not  available. 

Section  10  outlines  the  methods  for  the  transfer  and  Molniya-type  orbits. 

2  SATELLITE  DECAY  AS  A  'THREAT' 

There  is  a  widespread  perception,  among  the  public  and  particularly  the 
media,  that  the  decay  of  a  large  spacecraft  poses  a  threat  to  human  well-being. 
Consequently,  as  with  Skylab  1  in  1979,  there  can  be  a  voracious  demand  for 
forecasts  of  the  impact  point  or,  more  realistically,  of  the  track  over  the 
Earth  on  the  final  revolutions  in  orbit:  an  arcane  art  suddenly  becomes  of  public 
interest. 

The  'danger'  to  people  on  the  ground  from  pieces  of  satellites  falling  out 
of  the  sky  is  in  reality  minimal  by  comparison  with  the  other  hazards  of  life. 
There  are  about  two  satellites  per  week  from  which  pieces  survive  the  heating 
during  atmospheric  entry  and  tumble  to  the  ground,  and  the  impact  points  are 
scattered  almost  randomly  world-wide.  By  contrast  there  are  more  than  100 
incidents  per  week  when  aircraft,  or  parts  of  aircraft,  fall  out  of  the  sky,  and 
these  impacts  tend  to  be  concentrated  near  centres  of  population.  Also  some  of 
the  aircraft  are  very  large  and  some  are  filled  with  inflammable  fuel,  whereas 
the  satellite  fragments  are  inert  loops  (see  chapter  2  of  Ref  10  for  illus¬ 
trations).  So  far  as  is  known,  no  one  has  yet  been  injured  by  a  piece  of  a 
decaying  satellite,  and  the  probability  of  anyone  being  killed  in  the  next 
10  years  is  small,  whereas  more  than  50  people  per  year  on  the  ground  are  killed 
by  falling  aircraft  and  more  than  50000  per  year  by  traffic  accidents.  So  the 
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danger  of  satellites  to  people  on  the  ground  is  truly  negligible.  But  there  are 
no  limits  to  human  folly,  and  if  the  media  choose  to  label  a  particular  satellite 
as  'dangerous', ,  it  will  be  perceived  as  dangerous. 

Debris  from  satellites  failing  through  the  lower  atmosphere  does  pose  a 
greater  threat  to  aircraft  in  flight  m  that  region,  because  their  vulnerable 
area  is  so  large.  This  danger  is  still  very  small  by  normal  aircraft  safety 
standards,  but  it  may  be  piudent  to  divert  aircraft  known  to  be  flying  in  the 
vicinity  at  the  time  when  a  large  satellite  is  expected  to  decay. 

If  the  decaying  satellite  has  a  nuclear  reactor  as  a  power  source,  the 
scattering  of  radioactive  particles  in  the  atmosphere  would  have  some  long-term 
ill-effects.  However,,  the  satellites  are  small  and  the  ill-effects  are  likely 
to  be  ouch  less  than  those  caused  by  tests  of  atonic  bombs  in  the  atmosphere, 
by  accidents  in  nuclear  power  stations  and  by  the  small  amounts  of  radioactivity 
that  are  released  under  normal  operation  from  nuclear  installations. 

3  BASIC  METHOD  FOR  LIFETIME  ESTIMATION 

3. 1  Introduction 

There  are  two  possible  approaches  to  lifetime  estimation.  One  is  to  use  the 
observed  decay  rate  of  the  satellite  (with  an  approximate  value  of  eccentricity) 
to  assess  its  future  life.  The  other  is  to  calculate  the  lifetime  from  the 
current  orbital  elements  and  knowledge  of  the  size,  shape  and  mass  of  the 
satellite  and  the  current  air  density.  The  first  of  these  methods  is  much  to  be 
preferred,  because  the  second  method  requires  knowledge  of  the  mass  and  dimensions 
of  the  satellite,  and  these  are  not  accurately  known  for  about  90Z  of  the  satel¬ 
lites  launched.  Also  the  values  of  eccentricity  available  may  not  be  accurate 
enough  to  give  a  good  value  of  perigee  height;  and  even  when  the  perigee  height 
is  accurately  known,  the  current  density  at  perigee  is  tedious  to  calculate, 
because  it  depends  on  at  least  7  parameters;  and  even  then  will  only  be  accurate 
to  about  ±102.  In  this  section  the  first  method  is  covered,  but  an  outline  of 
the  second  method  is  given  in  section  9,  for  use  when  the  observed  decay  rate  is 
not  available. 

3.2  Practical  procedure 

The  most  widely  used  data  for  satellite  orbits  are  NORAD  2-line  elements 
issued  regularly  by  the  NASA  Goddard  Space  Flight  Center.  In  these  elements  the 
mean  motion  n  is  given  in  revolutions  per  day,  and  so  it  is  convenient  to 
express  the  lifetime  in  terms  of  n  ,  the  rate  of  change  of  n  expressed  in 
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revolutions/day  .  In  practice,  values  of  n  should  always  be  obtained  by  taking 
the  uifference  between  r jo  values  of  n  ,  and  dividing  b>  the  time  interval 
between  them,  which  needs  to  be  long  enough  to  ensure  a  change  of  at  least 
0.001  rev/day  m  n  .  The  value  of  h  given  m  the  2-line  element-  should  not 
be  used  (except  as  a  last  resort),  because  it  is  essentially  an  ’instantaneous’ 
value. 

The  basic  lifetime  L*  is  written  in  the  form 

L*  =  Q/n  ,  (1) 

where  Q  is  a  function  of  n  the  eccentricity  c  ,  the  atmospheric  density 
scale  height'  11  (which  depends  on  the  height  v  and  solar  activity),  and  the 
gradient  H  ,  that  is  dH/dy  which  we  denote  by  y  .  Values  of  H  and  y  for 
low  and  high  solar  activity  are  shown  m  Fig  !  The  values  are  from  the  COSPAR 

International  Reference  Atmosphere  1972 11  and  apply  when  the  solar  radiation 

-22  -2  -1 

energy  at  10.7,  cm  wavelength,  F ^  7  »  has  values  of  80  or  ISO  *  10  W  m  Hz 
(The  value  of  ^  is  usually  not  much  less  than  80  at  solar  nummum,  but  may 

greatly  exceed  150  at  a  strong  solar  maximum  ) 

The  variation  of  Q  with  e  for  high  and  low  solar  activity  and  for 
perigee  heights  of  200  or  600  km  is  shown  in  Fig  2.  For  e  *'  0  2  these  graphs 

are  based  on  equations  (5)  and  (6)  of  Ref  2,'  and  for  e  >  0.^  on  Ref  12:  there 

is  a  slight  discontinuity  at  e  =  0.2 

Although  solar  activity  does  not  appear  to  have  much  effect  in  Fig  2,' 
appearances  are  deceptive  because  lifetimes  most  often  have  to  be  calculated  for 
near-circular  orbits,  with  e  <  0.01,'  and  this  is  the  region  where  solar  activity 
has  most  effect.  More  detailed  diagram**  are  needed  for  accurate  calculation. 

The  first  of  these.  Fig  J,  shows  how  Q  varies  with  e  for  o  $  e  <  0.03 
for  various  values  of  n  ,  covering  perigee  heights  between  150  and  600  km.  lor 
perigee  heights  near  600  km,  the  change  from  low  to  high  solar  activity  increases 

Q  by  30Z  for  e  *=  0,  but  the  change  is  much  less  when  e  =  0.03  (about  5%). 

If  the  orbit  is  nearly  circular,  it  is  generally  better  to  work  with  the 
graph  of  Q  versus  n  for  0  $  e  <  0.01,  given  as  Fig  4 


+  H  is  a  measure  of  the  rate  of  decrease  of  density  p  with  height,  being  given  c 

by  H  -  -p/ (dp/dy) .  < 
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Fig  4  terminates  at  n  '  16  4  rev/d  ay  anil  does  nor  cover  the  very  lowest 
circular  orbits,  lor  a  circular  orbit  the  lifetime  is  given  by 


L* 


3Hs  n  ri 

2an 


(2) 


from  equations  (6.79)  and  (4.94)  of  Ref  j  In  this  equation  Hg  is  the  value 
of  H  at  height  Hg  below  perigee,,  n  is  given  by  equation  (4.88)  of  Ref  3 
and  a  is  as  usual  the  semi  major  axis  (nere  merely  the  radius  of  the  orbit). 

If  the  height  is  150  km,  winch  coriesponds  to  n  *  16.5  rev/day,,  Hg  -  13  km  for 
high  solar  activity  and  Hg  a  12  ka  for  low  solar  activiry,'  trom  tig  i.  also, 
if  the  lifetime  is  assumed  to  end  when  the  height  decreases  to  120  km, 

(ac  -  a^)  »  30  km  m  14.88),  and  B  =  0  1  ko  so  that  r,  ar0.95.  Equation  (2), 
then  gives  Q  *  0.047  for  high  solar  activity  and  Q  ®  0.043  for  low  solar  activity. 
These  values  are  indicated  by  circles  in  Fig  4. 

Solar  activity  has  most  effect  on  Q  for  near-circular  orbits  at  heights 
of  300-5C0  km,  as  is  obvious  from  the  wide  range  of  values  of  h  at  these  heights 

in  Fig  1.  For  example,  with  e  ■»  0  and  n  -  15.4  (corresponding  to  450  km  height). 

Fig  4  gives  Q  *  0.21  for  high  solar  activity  and  Q  *  0.16  for  low. 

When  0.03  $  e  $  0.2,  a  more  accurate  value  of  Q  can  be  obtained  by 
plotting  Q/e  against  e  for  various  values  of  n  ,,  as  shown  m  Fig  5.  For 

given  e  and  n  ,  the  variation  of  Q  with  solar  activity  never  exceeds  10%; 

and  the  variation  of  Q  with  n  ,  for  given  e  ,  is  also  nearly  always  less  than 
10%.  With  care.  Fig  12.4  can  be  read  with  an  accuracy  of  \Z,  which  is  slightly 
better  than  the  theory  on  which  it  is  based.  Although  Q  is  insensitive  to 
solar  activity  wnen  e  >  0.03,  the  lifetime  is  of  course  shorter  when  solar 
activity  is  high,  because  the  greater  atmospheric  density  leads  to  a  higher  value 
of  h  ,  and  hence  a  lower  value  of  Q/n  . 

Figs  3  to  5  are  the  same  as  the  corresponding  graphs  in  Ref  2.  For  e  >  0.2, 
12 

the  theory  has  recently  been  revised  ,  and  Fig  6  gives  the  variation  of  Q  with 
e  on  the  assumption  that  the  lifetime  is  reduced  by  a  factor  (1  -  0.1  u)  when 
y  is  taker  into  account,  this  is  close  to  the  result  obtained  in  Ref  12.  In 
using  Fig  6  it  should  be  remembered  that  when  e  >.  0.4,  quiu  large  lun  solar 
perturbations  iu  perigee  height  may  arise.  However,  if  -ne  amplitude  of  the 
lumsolar  perturbation  in  perigee  height  is  less  than  100  km,  and  if  tne  lifetime 
extends  to  several  cycles  of  perturbation,  L*  may  be  approximately  estimated  by 


calculating  the  average  value  of  n  over  a  cycle,  n  say*,  and  taking 
L*  *  Q/fi.  If  the  lunisolar  perturbations  are  large  and  not  easily  recognizable 
as  periodic,  numerical  integration  is  usually  needed  ,  though  analytical  methods 
are  still  sometimes  feasible  (see  section  !0). 

These  are  the  basic  techniques,  from  Ref  2  wich  some  additions.  We  now 
look  at  methods  for  taking  account  of  several  awkward  variations  that  sometimes 
spoil  the  accuracy  ot  the  basic  method. 

4  EFFECT  OF  OSCILLATIONS  IN  PERIGEE  HEICHT  CAUSED  BY  ODD  ZONAL  HARMONICS 

AND  ATMOSPHERIC  OBLATENESS 

Theory  for  these  two  effects,  which  both  produce  variations  in  perigee 
height  of  up  to  10  km,  has  been  developed  in  Ref  4.  If  L*  is  the  basic  life¬ 
time,  as  derived  in  section  3  from  spherical-atmosphere  theory,  the  true  life¬ 
time  L  may  be  expressed  as 

L  «=  L*  F(oz)F(ao)  ,  (3) 

where  F(oz)  is  the  factor  to  take  account  of  the  effect  of  odd  zonal  harmonics,, 
and  F(ao)  is  the  factor  that  allows  for  atmospheric  oblateness.  Expressions 
for  both  these  factors  are  derived  in  Ref  4. 

The  value  of  fi  in  equation  (1)  is  the  observed  current  or  (in  terms  of 
the  theory)  the  initial  value  of  n  ,  and  should  strictly  be  written  as  f»o  . 

For  simplicity,  the  suffix  o  has  been  omitted  m  section  3,  but  clarity  now 
demands  its  inclusion. 

The  gravitational  attraction  of  the  odd  zonal  harmonics  m  the  geopotential 
has  the  effect  of  producing  a  regular  variation  of  the  perigee  distance  r^ 
with  the  argument  of  perigee  u  (the  variation  is  virtually  sinusoidal  if 
e  >  0.01),,  The  amplitude  B  of  this  oscillation,  as  given  by  the  odd  zonal 
harmonics  of  Ref  13,  is  shown  in  Fig  7.  The  theory  of  Ref  4  requires  chat 
B/H  <0.3  and  shows  that  the  value  of  F(oz)  depends  on  the  parameter 


(*„>*) 


where  is  the  initial  (ie  current)  value  of  z  -  ae/H.  In  practice,  a  typical 

value  of  a/H  is  125,  so  that  z  ~  1  for  e  -  0.008  and  z  cr  iO  for  e  •  0.08.  The 


t  Averaging  1/A  is  not  recommended,  because  the  observed  value  of  n  can  be 
very  small  (or  even  slightly  negative  as  a  result  of  solar  radiation  pressure) 
when  perigee  height  is  at  its  maximum. 
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value  of  b  is  zero  for  zq  £  The  theory  of  Ref  4  is  not  valid  over  a  small 
range  of  inclinations  near  63.4  ,  where  B  becomes  very  large.  Since  B/H  must 
be  less  than  0.3,  the  maximum  permissible  value  of  B  depends  on  the  value  of 
H  :  generally  we  assume  H  >  33  km,  as  happens  if  the  perigee  height  exceeds 
180  km;  and  this  requires  that  B  <  10  km.  If  H  is  larger,  correspondingly 
larger  values  of  B  are  permissible. 

Fig  8  shows  how  F(oz)  varies  with  u>  when  b  -  0.2  for  various  values  of 

o 

u,  -  w  ,  where  <*>,  is  the  value  of  u  at  the  end  of  the  satellite's  life.  The 

L  O  L. 

thick  line  shows  the  value,  Fm  say,  when  u  goes,  truough  many  revolutions. 

The  thinner  lines  show  the  value,  F^  say,  when  u  compleces  only  a  fraction 
f  of  a  revolution,  so  that  l  -  u  l  -  2xf  .  If  w  goes  through  (N  ♦  f) 
revolutions,  F(oz)  may  be  taken  as  (NF^  ♦  fF^)/(N  ♦  f)  ,  or  may  be  calculated 
from  equation  (61)  of  Ref  4. 

It  is  also  shown  in  Ref  4  that  the  atmospheric  oblateness  factor  F(no) 

depends  on  and  X  :  here  X  depends  only  on  zq  and  is  shown  m  Fig  9, 

and  c'  is  the  initial  value  of  c*  «  cl«/l  ,  where  c  depends  on  the 
o  L  o 

ellipticity  of  the  atmosphere,  e  ,  and  is  given  by 


2H  po 


(5) 


with  i  as  the  orbital  inclination.  Also  I  denotes  the  Bessel  function  of 

n 

the  first  kind  and  imaginary  argument,  of  order  n  and  argument  z  ,  and  the 
variation  of  with  z  is  shown  in  Fig  10. 

Fig  11  shows  the  variation  of  F(ao)  with  when  c^  ■  0.2,  for  various 

values  of  u.  -  u  .  The  thick  line  shows  tb  results  for  u,  -  u  ■  ±Nr  ,  where 
L  o  L  o 

N  is  an  integer.  The  thinner  lines  show  the  results  when  u  does  not  complete 
a  half-revolution,  for  X  -  1  (corresponding  to  >  0.2):  the  broken  line  gives 
results  for  X  ■  0.8  (corresponding  to  zq  ^  8,  or  eQ  =?  0.06).  If  u  goes  through 
(N  +  f)  half  revolutions,  F(ao)  can  be  calculated  from  the  expression  given 
two  paragraphs  ago  for  F(oz)  ,  or  from  equation  (65)  of  Ref  4. 

The  values  of  b  and  c'  used  for  illustration  in  Figs  8  and  II,  although 

,  0  4- 

large,  are  within  the  possible  practical  range  for  high-inclination  orbits  ,  and 
the  Figures  show  that  changes  of  up  to  about  25%  in  lifetime  may  be  produced  by 
either  the  odd  zonal  harmonics  or  atmospheric  oblateness.  Both  correction  factors 
may  be  scaled  linearly  when  b  and  c^  differ  from  0.2. 

t  However,  can  only  be  as  large  as  0.2  if  c  »  0.3  and  zq  >  5. 


If  -  «o  is  sb.  11  enough,  F(oz)  and  F(ao)  differ  from  1  by  less 
than  10Z,  and  are  then  not  worth  evaluating:  this  happens  when  wl  "  w0  <  ^0° 
tor  F(oz)  or  when  <*>L  -  <  19°  for  F(ao)  . 

5  EFFECT  OF  THE  SEMI-ANNUAL  VARIATION  IN  AIR  DENS'.TY 

For  satellites  with  lifetimes  between  about  1  month  and  about  18  months, 
the  effect  of  the  semi-annual  variation  in  air  density  is  often  the  greatest 
source  of  error  ia  estimating  L  ,  Unfortunately  the  semi-annual  variation 
changes  appreciably  in  form  from  jear  to  year  and  cannot  at  present  be  accurately 
predicted.  The  only  practical  procedure  is  to  cake  a  standard  form  for  the 
variation,  and  to  hope  that  the  variation  which  is  in  progress  at  the  time  of 

the  prediction  is  not  too  seno  .sly  non-standard.  For  this  purpose  we  take  the 

14 

variation  shown  in  Fig  12,  which  is  an  average  over  several  years  and  applies 
at  a  height  of  aoout  250  km.  Most  satellites  with  lifetimes  of  1-18  months  have 
perigee  heights  between  200  and  300  km,  and  the  amplitude  of  the  semi-annual 
oscillation  does  not  vary  much  with  height  within  this  height  range,  so  the  use 
of  Fig  12  is  not  unreasonable.  However,  there  is  no  guarantee  of  its  accuracy 
in  future  yeais.  During  some  years,  errors  up  to  10Z  ’'r.  L  may  occur,  particu¬ 
larly  for  lifetimes  of  between  2  and  6  months;  during  other  years  the  variation 
may  be  near-standard,  and  errors  from  this  source  would  be  much  smaller.  It 
should  also  be  noted  that  the  semi-annual  variation  appears  to  be  greater  in 
the  northern  hemisphere  than  the  southern*^,  and  this  will  be  relevant  if  the 
satellite’s  perigee  is  expected  to  remain  in  the  same  hemisphere  until  decay. 

Again  we  define  L*  as  the  lifetime  in  the  absence  of  perturbations,  that 
is,  with  density  independent  of  time,  and  denote  the  ’semi-annual  correction 
factor’  by  F(sa)  ,  so  that  the  ’true’  lifetime  L  is  given  by 

L  -  L*F(sa)  .  (5) 

Then,  if  p  denotes  the  current  (initial)  value  of  the  ’standard  variation’  of 
o 

density  given  in  Fig  12,  and  if  p  denotes  the  mean  value  of  p  between  the 
initial  date  and  the  decay  date,  equation  (6)  of  Ref  4  shows  that 

F(sa)  »  pq/p  (7) 

Fig  13  shows  the  values  cf  F(sa)  for  lifetimes  up  to  160  days  and  a 
succession  of  current  dates  through  the  year,  on  the  assumption  that  the  semi¬ 
annual  variation  has  the  form  shown  in  Fig  12.  Usually  F(sa)  can  be  read 
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directly  from  Fig  13,  but  an  iterative  procedure  may  be  needed  if  F(sa)  departs 
greatly  from  1.0.  To  take  an  extreme  example,  if  the  basic  lifetime  is  calcu¬ 
lated  from  n  at  10  April,  and  L*  *=  80  days,  then  Fig  3  gives  F(sa)  =  1.25,, 
so  that  L  -  1.25L*  *  100  days.  The  appropriate  value  of  F(sa)  for  a  100-day 
lifetime  from  10  April  is  F(sa)  ■  1.31,  giving  a»  a  'second  iteration* 

L  -  I .31L*  -  105  days. 

An  alternative  plotting  of  F(sa)  against  the  time  of  year  for  a  specified 
lifetime  is  given  as  Fig  14  which  is  needed  for  lifetimes  greater  than  160  days. 
The  curve  tor  M  year'  applies  equally  well  lor  2,  3,  4,  ...  years. 

6  EFFECT  OF  THE  SCLAR-CYCLE  VARIATION  IN  AIR  DENSITY 

It  is  quite  useful  to  look  on  the  11-year  solar-cycle  variation  in  air 
density  as  a  slow-motion  version  of  the  semi-annual  oscillation.  Pursuing  this 
rather  loose  analogy,  we  see  that  the  basic  lifetime  L*  should  be  multiplied 
by  a  factor  F(sc)  ,  which  is  given  by 

F(sc>  -  cpl)/Jp  (8) 

where  p  is  the  current  (initial)  perigee  density,  and  p  the  mean  value  of 
po  '  p 

density  at  the  same  height  during  the  remaining  lifetime  L  . 

Unfortunately,'  solar  activity  is  ouch  more  irregular  than  the  semi-annual 
variation,  and  neither  the  timing  nor  the  amplitude  of  a  future  cycle  has  yet 
been  reliably  predicted  (except  by  luck),  though  there  are  several  confident 
predictions^  that  the  next  maximum  will  be  in  1990  or  1991,  with  a  maximum 
sunspot  number  of  1 1C  ±  10.  Fig  15  shows  the  sunspot  index  from  1610  to  1986, 
based  largely  on  Refs  19  and  °0.  The  irregularity  is  obvious,  and  the  possi¬ 
bility  that  solar  activity  might  return  to  the  very  low  levels  prevalent  m  the 
1 7th  century  cannot  be  ruled  out.  If  this  happened,  it  would  invalidate  life¬ 
time  predictions  made  on  the  assumption  that  the  mean  solar  activity  in  future 
cycles  will  be  similar  to  that  in  the  current  cycle  or  recent  past  cycles.  The 
lifetimes  would  then  be  much  longer  than  predicted. 

Until  the  solar  cycle  can  be  reliably  predicted,  it  is  onlj  possible  to 
adopt  empirical  and  imprecise  models.  If  the  satellite  is  expected  to  remain 
in  orbit  for  several  solar  cycles,  and  has  already  survived  one,  the  obvious 
procedure  is  to  use  the  increase  in  n  fin  say,  during  the  time  interval  fit 
between  one  solar  minimum  and  the  next  (or  between  one  maximum  and  the  next),  and 
to  take  L  »  Qfit/fin.  This  method  should  give  a  good  approximation  if  L  exceeds 
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about  three  cycles  and  if  subsequent  cycles  are  of  similar  strength  to  the  cycle 
used  in  calculating  fin  and  fit  .  If  the  subsequent  levels  of  activity  are 
higher  than  in  the  cycle  used,  the  real  lifetime  will  be  shorter  than  that  pre¬ 
dicted;  and  vice-ver3a. 

It  is  important  to  remember  that  decay  is  much  more  likely  to  occur  m  a 
year  when  solar  activity  is  high.  Indeed  solar  maximum  can  be  regarded  as 
'harvesr  time*  for  long-lived  satellites:  that  grim  but  sun-tanned  reaper  known 
as  *  solar  Max'  drags  them  down  in  hundreds.  Nearly  all  the  long-lived  satellites 
that  escape  the  cull  will  last  until  at  least  the  next  solar  maximum.  If  the 
predicted  lifetime  is  calculated  as  between  10  and  25  years,  we  need  to  adjust 
the  prediction  to  take  account  of  the  work  of  solar  Max.  For  example,  a  satel* 
lite  launched  at  solar  minimum  is  unlikely  to  have  a  lifetime  of  10  or  20  years, 
because  that  would  probably  imply  decay  near  solar  minimum,  which  is  unlikely. 

Ihe  lifetime  is  much  more  likely  to  be  near  5,  15  or  25  years,  and  one  of  these 
values  should  be  chosen  as  the  best  round-number  estimate  for  lifetime,  in 
preference  to  10  or  20  years. 

If  a  satellite  is  going  to  decay  within  the  current  solar  cycle,  the 
problem  is  even  more  difucult,  because  the  prediction  of  solar  activity,  and 
hence  of  air  density K  is  so  unreliable.  Only  a  very  rough-and-ready  method  can 
be  offered,  as  follows.  In  attempting  to  model  the  solar  activity  index  ^ 
over  a  solar  cycle,  a  sinusoidal  approximation  may  be  best;  but  the  density  has 
a  much  stronger  variation  than  Fjq  ^  (by  a  factor  of  up  to  30  at  500  km  height). 
If  is  the  density  at  perigee  height  at  solar  minimum,  and  the  value 

at  the  same  height  at  solar  maximum,  the  best  simple  analytical  representation 
of  the  density  at  this  fixed  height  seems  to  be 

P„  ■  0  +  (p„  -  »„)  sin4  (»t/P)  ,  (9) 

p  o  Me 

where  P  is  the  period  of  the  solar  cycle  and  t  is  the  ciue  measured  from 

solar  minimum.  In  practice  vhe  date  of  'solar  minimum'  is  ill-defined,  and  t 

should  be  measured  from  a  date  half-way  between  successive  solar  maxima:  the 

important  point  :s  to  try  to  ensure  that  the  maximum  density  p^  occurs  at  the 

right  time,  ie  that  t  *  JP  is  the  time  of  solar  maximum.  Equation  (9)  can  only 

be  called  a  desperate  expedient  rather  than  a  realistic  codel.  In  reality,  two 

successive  maxima  rarely  have  the  same  maximum  density  p„  ;  the  variation  is 

n 

asymmetrical,  the  rise  being  steepe  than  the  decline;  the  period  P  usually 
varies  between  one  cycle  and  the  next;  and  the  zero-date  for  t  is  subject  to 
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error  (for  rising  solar  activity)  because  the  date  of  the  next  maximum  may  be 
wrongly  forecast. 


If,  after  all  these  cautions,  we  still  adopt  equation  (9)  for 


value  of 


found  by  integrating 


3om)- 


p  between  time 
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In  equation  (10)  the  term  m  curly  brackets  is  in  practice  often  negligible, 
though  it  needs  to  be  retained  for  consistency  (to  ensure  that  p^  -*  p^o  as 
L-0).  P 

From  equations  (8),(  (9)  and  (10),,  the  explicit  equation  for  F($c)  is 


K5  •  V'.)  — rfr— 


•  <CM/Pft  -  1> 
2»(t  ♦  L) 


C'o'P) 


-KTr 


.  (ID 

When  large,  iteration  for  L  may  be  necessary  in  applying 

equation  (11).  Suppose,  for  example,  that  pM/pa  *  4  (appropriate  for  a  height 
of  about  330  km),  and  that  a  lifetime  estimate  made  at  solar  maximum  (t^  •  |P) 
gives  L*  ■  2  years.  Obviously  the  lifetime  will  be  longer  than  2  years,  and  if 
we  guess  L/P  -  1  for  simplicity,  equation  (11)  gives  F(sc)  *  1.30  or 
L  «  2.60  years.  If  P  •  II  years,  this  gives  L/P  -  0.236,  whence  F(sc)  ■  1.27 
from  (II)  and  L  *  2.54  years.  A  third  iteration  gives  L  «•  2  51  years. 

Fig  16  shows  the  variation  of  p^/p^  with  height  for  an  ’average*  solar 
maximum  (Fjq  ^  »  150,  giving  nr  1000  K  in  CTRA  1972 ),  and  also  for  a  high  solar 
maximum  (Fjq  ^  rs  250,  2  1300  K)  similar  to  that  of  1981.  At  solar  minimum, 

Tw  is  taken  as  700  K  for  the  purposes  of  Fig  16.  Most  satellites  with  lifetimes 
greater  than  half  a  solar  cycle  have  perigee  heights  greater  than  350  km,  so 
Fig  16  shows  that  in  practice  p^/pffl  usually  exceeds  4. 

Fig  17  shows  how  v(sc)  varies  with  t^/P  when  pM/pa  *  4,  for  selected 

values  of  lifetime  L  ,  aid  Fig  18  is  a  similar  diagram  far  pu/p  •  8.  As 

n  m 
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expected,  the  values  of  F(sc)  are  largest  when  the  lifetime  estimate  is  made 
at  or  a  little  after  sunspot  maximum  (t^  »  JP),  and  smallest  when  it  is  trade  at 
or  a  little  after  sunspot  minimum  (tQ  B  0) .  In  Fig  18,  F(sc)  can  be  as  large 
as  3  or  as  small  as  0.2. 

Values  of  up  to  20  are  possible  for  the  density  ratio,  as  Fig  16  shows,  when 

perigee  is  near  500  km.  So  the  variation  of  F(s»c)  with  p^/p  for  values  up  to 

20  is  shown  m  Fig  19,  for  L/P  •  1,2,3  ...  and  various  values  of  t  /P  . 

o 

Since  F($c)  is  so  strongly  dependent  on  the  density  ratio,  the  lifetime 
and  the  initial  poirt,  it  is  not  always  easy  to  interpolate  in  Figs  17  to  19, 
and  it  is  best  to  calculate  F(sc)  directly  from  equation  (11).  The  term  m 
square  brackets  in  equation  (11)  is  aero  if  L/P  is  an  integer,  so  that 

F(sc)  •  5  ♦"  HeJtJ  I1  4  (r  '  ')  K/P)|  (L/P  ■  l,:,3  ...)  . 

.  (12) 

This  form  may  also  be  used  as  an  approximation  whenever  L/P  >  10,  with  error 
less  than  ••6(pm/pq  "  1)2. 

In  Figs  17  and  18  the  value  of  F(sc)  often  departs  from  1.0  by  a  large 
factor:  so  lifetime  estimates  can  be  seriously  in  error  if  solar  activity  is 

ignored  or  wrongly  predicted.  The  error  is  greatest  for  lifetime  estimates  made 
21 

at  solar  minimum:  a  good  example  is  Skylab  1,  launched  in  1973,  for  which  NASA 
predicted  ».i  1974  a  total  lifetime  of  10  years,  whereas  UK  predictions  taking 
account  of  solar  activity  (correctly)  gave  6  years. 

EFFECT  OF  THE  DAY-TO-NIGHT  VARIATION  IN  AIR  DENSITY 
1  The  variation  of  density 

At  heights  above  200  km  the  air  density  is  highest  at  about  14  h  local 
time  and  lowest  at  about  04  h  local  time11.  The  variation  depends  on  latitude 
as  well  as  longitude,  however,  and  is  best  regarded  as  a  ’daytime  bulge*  m  the 
contours  of  density,  having  its  maximum  near  the  equator  (or,  more  accurately, 
near  the  sub-solar  latitude)  and  lagging  behind  the  sub-solar  longitude  by  about 
30°.  Consequently  the  perigee  of  a  satellite  of  high  inclination,  as  it  slowly 
travels  round  the  orbit,  sometimes  avoids  the  daytime  bulge  by  'slipping  past' 
at  high  latitude,  and  sometimes  goes  through  the  centre  of  the  bulge.  The  track 
depends  on  the  right  ascension  of  the  node  and  the  time  of  year.  The  result  is 
that  the  day-to-night  variation  in  perigee  density  experienced  by  the  orbit. 
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although  nearly  constant  in  period,  varies  in  amplitude  as  the  right  ascension 
of  the  node  steadily  precesses. 

The  density  at  perigee  depends  primarily  on  the  geocentric  bulge-perigee 
angle  as  specified  in  Ref  22.  For  orbits  of  high  inclination,  the  variation 

of  6  ditfers  free  cycle  to  cycle  (as  seen  in  Figs  5  to  8  of  Ref  14,  for 
P  o 

example),  and  a  regw* .  variation  of  can  only  be  assumed  if  i  <  35  .  How¬ 
ever,  this  simplest  of  variations  doer  also  serve  as  a  fair  approximation,  on 
average ,  for  satellites  of  high  inclination,  if  they  experience  enough  day-to- 
night  cycles,  although  the  average  amplitude  is  then  generally  about  0.6  of  the 
maximum  possible,  because  the  perigee  rarely  samples  the  full  maximum  and  more 
often  crosses  the  bulge  on  its  lover  slopes.  So  a  correction  factor  based  on 
the  assumption  that  the  variation  is  regular  will  be  useful  for  satellites  of 
low  inclination  and  for  long-lived  satellites  of  higher  inclination. 

The  most  realistic  simple  approximation  for  the  variation  of  density  p 
with  angular  distance  $  from  the  centre  of  the  daytime  bulge  appear  (from 
section  3.2  of  Ref  22)  to  be 

P  -  pm  +  (PM  "  Pm>  cos  1*  » 

where  pm  and  p  are  now  the  maximum  and  minimum  densities  in  the  day-to-night 
variation.  Applying  the  above  equation  for  perigee  (with  p  ■  p^  and  >  ■  $p) , 
and  assuming  ♦  varies  linearly  with  time  t  measured  from  the  minimum-density 

.  t  ? 

point  ,  we  have 

Qn  -  P  +  (PM  -  P  )  sin4  (at/P)  ,  (13) 

p  m  n  m 

where  the  period  P  is  now  that  of  the  day-to-night  cycle  rather  than  the  solar 
cycle.  Equation  (13)  is  the  same  as  equatior  (9);  consequently  p  is  given  by 
equation  (10),  and  F(dn)  by  equation  (II). 

23 

The  period  P  of  the  day-to-night  cycle  may  be  written  as 

P  -  - - 3^2 -  days  ,  (14) 

4.982(R/a) (Scos  i  -  2cos  i  -  1)  -  0.986 


t  Formally,  this  implies  that  fy,  •  *(l  -  2t/P)  .  In  reality,  ♦p  reverts  to 
positive  values  when  0.5  <  t/P  <  I .0,  but  its  sign  is  of  no  consequence  in 
evaluating  cos4(J$  )  . 
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if  (1  -  e  )  is  taken  as  1-  Values  of  P  given  by  equation  (14)  art  plotted 
against  i  for  R/a  ■*  0.8  and  0.9  m  Fig  20.  This  diagram  is  difficult  to  read 
when  the  period  exceeds  1  year  and  values  of  are  given  as  an  alternative. 

If  P  is  in  yeais  and  in  degrees/day,  P  *  0.986/$^. 

The  appropriate  amplitude  foi  the  ratio  of  maximum  to  minimum  density, 
p^/p^  ,  in  the  day-to-night  variation  indicated  m  Fig  2 1 wnere  the  values 
of  density  are  taken  from  C1RA  1972 
for  for  'low  solar  activity',  a 

for  for  'high  solar  activity' . 

In  practice,  satellites  with  perigee  heights  greater  than  about  500  km 
usually  have  lifetimes  which  extend  over  many  cycles  of  the  day-to-nigh‘_  variation, 
and  it  is  often  simplest  then  to  calculate  n  as  an  average  over  one  cycle  and 
not  to  use  F(dn)  at  all.  If  the  F(dn)  method  is  used  mainly  for  perigee 


most  likely  to  arise  in  practice. 


By  deriving  equation  (13),' ve  have  tacitly  implied  that  the  density  at 
perigee  will  be  a  crucial  parameter.  This  is  likely  to  be  true  if  the 
eccentricity  is  high  enough,  but  not  for  a  circular  orbit.  Thus  we  need  to 
divide  orbits  on  the  basis  of  their  eccentricity and  e  *  0.03  proves  a  con¬ 
venient  dividing  line. 

7.2  Orbits  with  c  >  0.03 

It  is  shown  in  the  Appendix  that,(  for  orbits  with  e  >'0.03  moving  in  an 
atmosphere  having  the  simple  sinusoidal  variation  of  p  with  $  specified  in 
Ref  22,  the  factor  F(dn)  that  takes  account  of  the  day-to-night  variation  m 
air  density  may  be  written  as 


F(d,.)  .  %o/pp  (15) 

where,  as  before,  p  is  the  initial  (current)  value  of  p  and  p  is  the 
po  P  P 

mean  value  of  density  at  the  initial  perigee  height.  Equation  (15)  will  be 

taken  as  the  'standard'  Hdn)  ,  and  it  is  of  course  the  same  as  equation  (8)  for 
F(sc)  .  Equation  (15)  does  have  limitations,  however:  its  accuracy  would  be 
dubious  if  the  lifetime  was  less  than  one  cycle  of  ,  because  the  model  is 
then  probably  inadequate.  Also  it  applies  strictly  only  for  orbits  of  incli¬ 
nation  less  than  35°,  or  for  higher-mclination  orbits  over  several  cycles  of 

$  ;  and  of  course  e  must  exceed  0.03. 

P 
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We  assume  that  equation  (15)  still  holds  when  p  is  represented  by  a  aore 
accurate  moael,  such  as  equation  (13),  Then,  within  the  provisos  of  the  previous 
paragraph,  the  variation  of  F(sc)  with  tQ/P  in  Fig  17  applies  equally  veil  to 
F(dn)  ,  although  the  broken  curves  may  be  unreliable  because  they  are  for  values 
of  L/P  less  than  I . 

In  using  Fig  17  to  obtain  F(dn)  ,  the  appropriate  period  P  is  of  course 
that  of  the  day-to-night  cycle,  and  tQ/P  represents  the  fraction  of  that  cycle 
which  has  been  completed  at  the  date  when  the  lifetime  estimate  is  made.  Thus 
if  t^  »  0,  perigee  is  at  the  minimum-density  point;  if  tQ  »  lP,  perigee  is  moving 
towards  the  maximum  density  and  is  half-way  there;  if  »  JP,  perigee  is  at  the 
maximum-density  point;  tf  tQ  •  $P,  perigee  is  half-way  on  its  journey  from 
maximum  to  minimum  density. 

If  p„/p  '<  4.  as  is  likely  m  practice,  the  values  of  F(dn)  obtained  from 
n  m 

Fig  17  can  be  scaled  down  by  use  of  Fig  IS.  A  more  accurate  procedure  is  to 

calculate  F(dn)  from  equation  (12),  if  L/P  •  1,2,3  ....  If  p^/p^  <  2  K  a 

linear  variation  of  F(dn)  with  p^/o^  may  be  assumed.  Thus  if  F(dn,2) 

denotes  the  value  for  Pu/p  *  2,  the  value  of  F(dn)  for  any  other  p,./p_  less 
n  m  Mm 

than  2  may  be  taken  as 

F(dn)  -  1  »  |oH/Pn  -  l||F(Jn,2)  -  l|  <'  «pm/pb«2)  •  (|6> 

7.3  Orbits  with  e  <  0.03  (z  <  3) 

Values  <.f  F(dn)  for  e  <  0.03  are  derived  m  the  Appendix,  where  it  is 

shown  that  the  standard  value  of  F(dn)  should  be  scaled  down  towards  1  by  a 

factor  I,(z  )/ln( 2  )  ,  which  decreases  from  0.81  at  z  -  3  to  0  at  z  »  0  (see 
I  o  u  o  o  o 

Table  Al).  Thus 

{FCdn)}^  -  '  *w{FWn>  '  <l7> 

If  pM/p  <  3,  as  is  likely  in  practice  for  most  high-inclination  orbits.  Fig  19 
shows  that  0.6  <  F(dn)  <  1.7.  The  factor  |F(dn)|2  may  therefore  still  be 
significant  down  to  quite  low  values  of  e  :  to  keep  it  within  10Z  of  1 .0 
requires  z  <  0.3,  that  is,  an  eccentricity  below  about  0.003. 

7.4  Review  of  results 

It  should  be  emphasized  that  the  day-to-night  variation  is  a  very  difficult 
effect  to  treat,  and  the  results  have  many  limitations  and  approximations. 
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Strictly,  the  results  for  F(dn)  are  realistic  only  for  i  '<  35°  or  for  satellites 

of  higher  inclination  with  lifetimes  extending  over  several  day-to-night  cycles. 

For  the  latter  the  maximum  density  is  rarely  sampled  and  a  'scaling-down  factor' 

of  about  C.6  is  recommended  to  find  the  effective  value  (PU?P  )  ee  for  use 

n  Q  etr 

in  Fig  17  or  19.  Thus 


Wef i  •  1  •  <1S> 

where  on  the  right-hand  side  is  found  from  Pig  21.  After  this  correction 

has  been  made  (when  necessary),  the  values  of  F(dn)  can  be  read  from  Fig  17  or 
19,  and  then  modified  using  equation  (17)  if  zq  <  3  (eft  <  0.03). 

Finally,  it  should  be  noted  that  if  L/P  “■  1,2,3  ...  ,  or  if  L  >  10P,_  the 
value  of  F(dn)  may  be  tound  from  equation  (12),  with  the  reduced  value 
l£  1  A5°- 

8  REVIEW  OF  THE  SOURCES  OF  VARIATION  IN  LIFETIME 

8.1  Sumnary  of  procedure  for  applying  the  five  correction  factors 

Fro*a  the  results  of  sections  4  to  7,  the  full  procedure  for  predicting 
lifetime  is,  first,  to  calculate  the  basic  lifetime  L*  using  Figs  3  to  6, 
and  then  to  apply  up  to  five  correction  factors,  the  true  lifetime  L  being 
expressed  as 


L  «  L*  F (02)  F(ao)  F(sa)  F(sc)  F(dn)  .  (19) 

The  odd  zonal  correction  factor,  F(oz)  ,  is  given  by  equation  (61)  of  Ref  4,  and 
here  by  Fig  8,  with  Figs  1  and  7  as  auxiliaries.  The  second  correction  factor,, 
F(ao)  ,  takes  account  of  atmospheric  ob lateness,  being  given  by  equation  (65)  of 
Ref  4,  and  being  shown  here  m  Fig  11,  with  Figs  9  and  10  as  auxilianes.  The 
third  correction  factor,  F(sa)  ,  takes  account  of  the  semi-annual  variation  in 
air  density  and  is  shown  in  Figs  13  and  14.  The  solar-cycle  factor,  F(sc)  , 
is  defined  by  equations  (8)  to  (1!)  and  plotted  in  Figs  17  to  19,  with  Fig  16 
as  an  auxiliary.  The  final  (and  least  satisfactory)  correction  factor,  F(dn)  , 
provides  an  indication  of  the  effects  of  the  day-tt -night  variation  in  air 
density.  Standard  values  of  F(dn)  ,  which  are  the  same  as  for  F(sc)  ,  are 
given  in  Figs  17  and  19,  though  now  with  Figs  20  and  21  as  auxiliaries:  strictly, 
the  standard  F(dn)  applies  only  when  (a)  perigee  goes  through  at  least  one  day- 
to-night  cycle,  (b)  the  inclination  is  less  than  35°,  and  (c)  the  eccentricity 
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is  greater  than  0.03.  The  correction  is  smaller  for  orbits  of  high  inclination 
and  for  eccentricities  less  than  0.03,  as  explained  m  section  7. 

It  is  rarely  that  the  five  factors  are  significant  simultaneously,  but  all 
can  have  important  effects  on  the  lifetime.  If  we  take  any  change  greater  than 
10%  to  be  significant,  F(ao)  may  be  of  importance  for  lifetimes  greater  than 
2  weeks,  and  F(oz)  and  F(sa)  for  lifetimes  greater  than  1  month.  All  these 
three  factors  generally  have  maximum  effects  of  about  30%'.  The  effect  of 
F(sc)  can  be  significant  for  lifetimes  greater  than  1  year  and  is  very  large 
for  satellites  with  perigee  heights  near  500  km  having  lifetimes  of  5-20  years: 
Fig  18  shows  that  the  lifetime  estimate  may  have  to  be  reduced  by  a  factor  of 
up  to  5  if  made  at  solar  minimum,  or  increased  by  a  factor  of  up  to  3  if  made 
at  solar  maximum.  The  effect  of  F(dn)  is  small  if  the  perigee  height  is  below 
250  km,  but  can  become  large  (>50%)  for  satellites  having  lifetimes  of  several 
years  and  perigee  heights  above  250  km. 

8.2  Other  effects 

Although  five  important  variations  are  covered  in  the  previous  paragraph, 
the  effects  of  short-term  density  variations  due  mainly  to  solar  activity  have 
not  been  brought  into  the  scheme  of  calculation  because  they  are  as  yet  unpre¬ 
dictable.  A  useful  (though  arbitrary)  division  ot  these  is  into  four  categories: 

(1)  Sudden  changes  m  density  occurring  within  a  few  hours,'  attribu¬ 
table  to  transient  outbursts  of  solar  activity  and  closely  correlated 
with  geomagnetic  storms.  For  lifetimes  greater  than  I  week,  the  tiroes  of 
such  disturbances  should  be  avoided  when  calculating  n  .  If  the 
satellite's  lifetime  is  only  a  few  days,  it  is  wise  to  look  at  the  pre¬ 
dicted  values  of  the  geomagnetic  index  and  to  make  allowance  for  any 
unexpected  disturbances.  These  sudden  changes  linked  with  geomagnetic 
storms  are  usually  most  prevalent  in  the  year  or  two  after  solar  maximum. 

(ii)  Monthly  changes.  The  Sun  rotates  relative  to  the  Earth  once  every 
27.3  days,  and  sunspot  groups  that  persist  on  the  Sun's  face  for  several 
rotations  tend  to  produce  a  27-day  recurrence  in  upper-atmosphere  density 
variations.  If  such  recurrences  have  arisen  in  preceding  months,  it  is 
best  to  evaluate  n  over  about  27  days,  if  the  lifetime  exceeds  about 
6  months.  For  shorter  lifetimes,  fi  should  if  possible  be  evaluated  at  a 


t  The  effect  of  F(oz)  can  be  larger  at  inclinations  close  to  63. A0,  which  are 
excluded  from  the  theory. 
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level  of  solar  activity  that  is  likely  to  be  the  average  over  the  remaining 
lifetime  -  a  procedure  that  calls  for  some  inspiration  (or  guess-work) 

(lii)  Irregular  short-term  variations.  When  all  allowances  have  been 
made  for  the  effects  of  solar  activity,  day- to-night  variations,  etc,  there 
usually  remains  an  irregular  variation  m  density,,  and  hence  n  having  an 
amplitude  of  about  *102  and  a  timescale  r arfging  between  a  day  and  a  week; 
and  this  produces  errors  of  up  to  102  in  lifetime  estimation  because  of 
the  variability  in  i.  . 

(iv)  The  possibility  of  large  changes  m  solar  activity.  In  the  two 
or  three  years  before  and  after  a  strong  solar  maximum,  solar  activity  may 
undergo  quite  large  and  unexpected  variations  which  can  cause  larger  errors 
in  lifetime  prediction,  particularly  for  lifetimes  between  about  6  months 
and  2  years.  These  errors  are  less  likely  to  arise  in  the  2  years  before 
and  after  solar  minimum. 

When  the  satellite  is  within  a  few  revolutions  of  decay,  two  further 
sources  of  variation  are  liable  to  cause  errors  in  lifetime  estimation. 

The  first,  which  applies  to  satellites  of  all  shapes  and  si?es,;  arises 
because  the  drag  coefficient  does  not  remain  constant  as  assumed  in  the  theory: 
instead  there  is  a  decrease  in  the  drag  coefficient  as  the  perigee  descends  to 
heights  where  frec-molecule  flow  no  longer  prevails.  At  heights  of  150-400  km 

">4 

the  drag  coefficient  is  often  taken”  to  be  about  2.2,  but  the  value  begins  to 

decrease  when  perigee  descends  below  130  km  and  the  drag  coefficient  is  reduced 

25 

to  about  half  its  frec-molecule  value  by  the  time  perigee  has  descended  to 
90  km.  At  this  stage  there  has  already  been  some  extension  of  the  lifetime,  and 
the  subsequent  time  in  orbit  will  be  double  what  the  theory  predicts.  The  net 
effect  varies  between  one  satellite  and  another,  but  a  rough  rule  of  thumb  is  to 
add  one  revolution  to  the  lifetime  for  predictions  made  in  the  last  24  hours  of 
the  satellite’s  life.  This  very  rough  rule  is  most  appropriate  for  quite  large 
satellites  with  mass;area  of  order  100  kg/o  and  having  orbits  of  low  eccen¬ 
tricity.  The  effect  is  smaller  if 

(a)  the  satellite  is  small  (because  it  remains  in  free -molecule  flow 
down  to  a  lower  altitude), 

(b)  the  oass/area  is  lower  (because  the  remaining  life  is  then  shorter),  or 

(c)  the  orbit  is  of  high  eccentricity  (because  perigee  may  already  be  below 
HO  km  some  days  before  decay). 

Fig  26  of  Ref  2  covers  this  effect  graphically. 
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The  second  source  of  error  applies  to  all  non-sphencai  satellites  and 
arises  from  changes  in  attitude,  which  alter  the  drag  or  create  lift.  The 
simplest  variation  is  a  change  in  the  mode  of  rotation  of  the  satellite,  produced 
by  aerodynamic  forces  at  perigee:  at  worst,  a  tumble  can  be  converted  to  arrow¬ 
like  flight,  thus  considerably  reducing  drag  and  extending  the  lifetime;  tms  is 
unusual,  but  changes  of  ±20%  in  drag  probably  occur  quite  often,  especially  for 
satellites  of  complex  shape  (such  as  those  with  solar  panels).  A  more  drastic 
variation  in  the  mass/area  occurs  if  part  of  a  satellite  is  torn  off  by  the 
effects  of  aerodynamic  heating  and  forces  at  perigee:  the  remaining  structure 
then  in  effect  becomes  a  new  satellite,  needing  a  new  lifetime  estimate,  based 
on  subsequent  orbital  data,  which  will  not  exist  unless  the  satellite  survives 
for  several  more  revolutions.  A  third  problem,  even  more  intractable,  arises 
when  the  attitude  of  the  satellite  as  it  approaches  perigee  is  such  that  it 
generates  upward  lift,  tends  to  'skip  off*  the  lower  atmosphere  and  survives  for 
a  few  more  revolutions.  (If  the  lift  is  negative,  it  may  'dive  in'  earlier  than 
predicted.)  It  is  easy  enough  to  think  of  other  scenarios  disastrous  to  math¬ 
ematical  prediction,  for  example,  aerodynamic  heating  at  perigee  may  be 
sufficient  to  provoke  an  explosion  in  residual  propellants  remaining  in  a  rocket. 

9  LIFETIME  IN  ~ RMS  OF  p  AND  6 
9. 1  Introduction 

The  basic  method  of  section  3  requires  a  value  for  the  current ly-cb served 
decay  rate,  A  .  Sometimes  this  will  not  be  available,  but  it  may  still  be 
possible  to  calculate  the  satellite's  lifetime  if  the  mass  and  dimensions,  and 
the  perigee  height,  are  quite  accurately  known.  In  the  calculations  a  spherical 
atmosphere  will  be  assumed:  thus  it  is  L*  that  will  be  found,  and  th-  correc¬ 
tions  F(oz)  etc  can  then  be  applied  as  in  section  8. 

The  lifetime  will  be  expressed  in  terms  of  the  density  at  perigee,  , 
and  the  arca/mass  parameter,  6  ,  defined  by 

5  -  FSCp/m  .  (20) 

Here  m  is  the  mass  of  the  satellite,  S  its  mean  cross-sectional  area  perpen¬ 
dicular  to  the  direction  of  motion,  C_  is  the  drag  coefficient,  for  which  a 
24  u 

value  of  2.2  is  recommended  for  perigee  heights  of  150-400  km,  and  F  is  a 
factor  to  allow  for  atmospheric  rotation,  which  is  usually  between  0.9  and  1.1 
and  is  shown  in  F*g  2.5  of  Ref  3. 
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9.2  Orbits  with  e  <  0.2 

With  a  spherical ly-sytmmitncal  attmisphcre,  equation  (4.22)  ot  Rc£  3.  v.th 
a  .  a  ,  shows  that  we  may  write  the  change  ha  in  a  during  one  revolution  as 
o 

ha  "  '  2»aInp«  exp(-t>{l0  *  2el,  ♦  0(e2)}  ,  <“  ‘  °'2>  W 


after  dropping  the  suffn  o  from  a  ,  Pp  and  a  . 
by  Kepler's  third  law,  & n  -  -3-iha/ l2a)  ;  and,;  as  it 
period  T  wc  have 


2  3 

Since  n  a  *  constant 
is,  equal  to  the  orbital 


An  3nAa  m  _  3&a  (22) 

AT  2aAt  2ar 


tf  n  is  expressed  in  revolutions  per  day  (rather  than  rad/day),,  so  that 
n  „  |/T.  From  equations  0).  (22)  and  <20,  we  have 


L* 


X(e,H) 


(23) 


vhete 


X(e,H) 


QT  exp  2 

10(S)  ♦  delj (z)  *  <He2) 


(2M 


for  e  '<  0.2.  Values  of  X  are  plotted  against  e  for  two  values  of  H/a  in 
fig  22.  If  a  ,  e  ,  T  ,  P  and  11  are  known,  equation  (23)  gives  the  life¬ 
time  for  e  -<  0.2.  Approximate  values  of  Pp  for  two  levels  of  solar  activity 
are  available  from  Fig  23,  but  a  more  detailed  calculation  using  CiM  « 

advisable.  The  same  advice  applies  for  H  ,  for  which  F,g  .  gives  approximate 

values . 

It  is  often  useful  to  eliminate  T  from  equation  (23)  by  writing^ 

I  -  2.(a3/:  ’)!  ,-  where  v'  is  the  gravitational  constant  (338bOO  km  s  for 
the  Earth).  Equation  (23)  then  becomes 


L* 


X(c,H) 


(">5) 


The  function  X(e,H)  depends  on  c  and  H/a 
f  equation  (4.56)  of  Bof  3,  QT  *  ef(e 


only,  because,  in  the  notation 
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When  e  =  0,  Che  lifetime  can  be  obtained  from  equation  (4.89)  of  Ref  3, 
which,  on  dropping  the  suffix  and  writing  T  *  2*(a^/y,)J,  gives 

L*  «  - !_  .  (26) 

(y  *a) *p6 

For  satellite  heights  of  200-IC00  km,  (y’a)^  varies  by  less  than  32  from  its 

15  2  -1 

value  at  600  km  height,  which  is  4.55  x  10  m  day  .  So  equation  (26)  gives 
the  lifetime  in  a  circular  orbit  3S 

„  0.22H  .  .... 

L*  ~  — re —  days  (27) 

!G  3p6 

3  2 

where  H  is  in  metres,  p  in  kg/m  ,  and  5  in  aT /kg  If  we  further  take 
0^  *  2.2  and  F  ■  1 so  that  6  *  2.2  (S/m),  equation  (27)  reduces  to 

L*  ~  (f)— 7 7T  days  '  (28) 

Vs'  10  8p 

~12  3  3 

When  o  *  10  kg/m  (1  ng/m  )  t/  which  from  Fig  23  is  the  value  at  about  550  km 
height  when  solar  activity  is  high,  and  H  »  70  km  (from  Fig  1),  equation  (28) 

2 

gives  L*  cs  700  days  if  m/S  -  100  kg/m  v  a  typical  value.  At  the  same  height 

.  “14  3 

with  low  solar  activity,  the  appropriate  o  would  be  about  8  x  10  kg/ta  , 

from  Fig  23,  and  H  would  be  about  56  km,  from  Fig  1:  so  L*  -  70  m/S,  or 
2 

20  years  if  m/S  *  100  kg/m  .  This  value  of  L*  is  unrealistic  because  low  solar 
activity  has  (m  the  past  250  years)'  never  persisted  for  20  years;  but  the  calcu¬ 
lation  does  show  that  lifetimes  in  circular  orbits  can  be  up  to  10  times  longer 
at  solar  minimum  tha  i  at  solar  maximum  for  heights  near  500  km. 

For  circular  orbits,  X(e,H)  -  1.5  H/a,  as  comparison  of  equations  (25)  and 

(26)  shows:  the  lifetime  is  longer  when  H  is  larger,  because  the  increase  of 

density  as  time  goes  on  is  less  rapid  when  11  is  larger.  For  e  *  0.01,  however, 

Fig  22  shows  that  Xte,H)  is  larger  when  H  is  smaller,  this  is  because  the 

density  above  perigee  decreases  more  strongly  when  H  is  smaller,  thus  giving 

rise  to  a  smaller  A  for  given  p  and  e  .  When  z  is  large,  X(e,H)  is 
P  l 

approximately  proportional  to  (a/H)*  for  given  e  .  Inese  variations  with  H 
are  of  course  relatively  small:  for  given  e  ,  it  is  mainly  0^6  that  controls 
the  lifetime  in  equation  (25). 
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9.3  Orbits  with  e  >'  Q.2 

The  results  so  far  have  been  for  e  s  0.2.  When  e  >  0.2,  equation  (2 1 )  has 
to  be  replaced  by  the  high-eccentricity  version,  equation  (4.107)  of  Ref  3,  which 
gives 

|  212 

4a  -  -  a2c  *  C>,  *  0(0. 003) i  (e  >  0.2)  .  (29) 

p  '*  '  (1  -  ey11  *  ' 

Then  L*  is  given  by  equation  (23)  with 


X(e,H) 


QT|2aa(l  -  u) | 


(1  +  a) 


372 


(e  >  0.2) 


(30) 


Writing  QT  -  eF(e)  in  the  notation  of  equation  (4.133)  cf  Ref  3,  and 

z(l  -  e)  »  r  e/H,  we  have 
P 


X(e.K) 


(e  >  0.2) 


(31) 


Using  the  values  of  F(e)  given  in  Fig  4.10  of  Ref  3,  and  taking  H/r^  s  0.01, 
we  find  X(e,H)  increases  from  i.6  at  e  *  0.2  to  4  4  at  e  *  0.4  and  13.2  at 
e  -  0.8, 


9.4  Results  for  average  solar  activity 


Fig  24  is  a  graph  of  satellite  lifetime  as  a  function  of  perigee  height 
and  e  ,  the  values  of  peri^ae  density  0^  and  scale  height  H  being  taken 

from  CIRA  1972  with  *  900  K,  which  is  approximately  the  average  over  a  solar 

°°  .  2 
cycle.  The  quantity  plotted  in  Fig  24  is  L*(S/m);  so,  if  m/S  «  100  kg/m  ,  the 

lifetime  m  years  is  100  times  the  value  shown.  In  drawing  Fig  24  (which  vs 

from  Ref  2),  it  was  assumed  that  CD  =  2.2  and  F  »  0.91.  If  the  lifetime  obtained 

from  Fig  24  is  less  than  10  years,  the  solar  activity  is  unlikely  to  be  'average*, 

and  the  mean  solar  activity  over  the  predicted  lifetime  ought  to  be  used.  For 

such  short  lifetimes  the  best  technique  is  to  calculate  L*  with  the  current 

perigee  density, (  from  equation  (23)  or  (27),  and  then  to  correct  for  solar 

activity  by  the  method  already  given  m  section  6.  The  other  corrections  to 

L*  can  be  applied  whenever  they  are  appropriate  and  appreciable 


i 


\ 
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1 0  HIGH-ECCENTRICITY  ORBITS  WITH  LUNISOLAR  PERTURBATIONS 

10.1  Introduction 

The  opthods  of  sections  3  to  9  are  valid  for  eccentricities  greater  than 
0.2,  if  lunisolar  perturbations  m  perigee  height  are  small  enough.  When 
e  >  0.4,  however,  lunisolar  perturbations  usually  have  a  significant  effect  and 
need  to  be  computed.  If  the  lunisolar  perturbations  produce  a  variation  which  is 
recognizably  periodic,  the  lifetime  can  be  estimated  by  using  a  value  of  n 
averaged  over  a  cvcle;  if  the  perturbation  in  perigee  height  is  large  (>50  km) 
and  not  recognizably  periodic,  it  is  usually  necessary  to  resort  to  numerical 
integration.  However,  it  is  worth  trying  to  devise  analytical  methods  for  some 
of  the  heavily-populated  “ypes  of  high-eccentricity  orbit,  and  two  such  methods 
will  now  be  described,  for  low-inclination  transfer  orbits  to  geosynchronous 
height  and  for  Molmya-type  orbits. 

10.2  Low-inclination  transfer  orbits  (i  »  25°,  eo=  0.73) 

I  >sl  of  the  ma.4y  geosynchronous  satellites  now  in  orbit  have  made  their 
way  to  the  geosynchronous  height  of  36000  km  via  a  transfer  orbit  with  perigee 
height  near  200  km  and  apogee  height  36000  km.  Such  transfer  orbits  have 
eccentricities  close  to  0.73,  and  most  have  orbital  inclinations  between  20° 
and  30°  A  discarded  rocket  usually  remains  in  the  transfer  orbit  and,  as 
there  are  now  more  than  a  hundred  ot  these,  an  approximate  analytical  method  for 
determining  their  lifetimes  has  been  developed  . 


In  the  absence  of  lunisolar  perturbations  the  lifetime  L*  is  given  by 
equations  (25)  and  (31),  and,  on  inserting  the  numerical  values  of  a  and  e 
for  a  transfer  orbit,  and  taking  CD  «  2.2  and  F  »  0.925,  we  have 


23.6 


(  ioos) 


years 


(32) 


.  2  .  3 

with  m  in  kg,  S  in  m  ,  o  m  kg/m  and  H  in  km.  The  variation  of  L* 

P  2 

with  perigee  height  for  a  standard  value  of  m/S  ■  100  kg/m  is  shown  in  Fig  25, 
2 

If  m/S  differs  from  100  kg/m  ,  the  lifetime  should  of  course  be  multiplied  by 
m/lOOS  . 

The  average  period  of  the  lunisolar  perturbation  in  perigee  height  for  a 
transfer  orbit  is  about  9  months,  though  the  amplitude  and  period  of  the 


t  The  exceptions  (i  =10°)  are  satellites  launched  by  Ariane. 
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perturbation  vary  appreciably  from  cycle  to  cycle.  If  the  satellite  survives 
several  cycles  of  the  lunisolar  perturbation,  we  can  calculate  the  perturbation 
approximately,  determine  the  mean  perigee  height  and  take  that  as  the  'perigee 
height'  in  Fig  25.  It  can  be  shown^  that  if  A  is  the  solar  longitude  and 
i  ■  25°.  the  arithmetic  mean  perigee  height  y^  is  give  »  by 

y  e  y  ♦  33  cos<2(ft  ♦  u  *  A  )>  +  26  cos  (ft  +  2w  )  “21  cos<2(fl  +  w  )>  Km 
7p  7po  (  o  o  o  /  o  o  (oof 


(33) 


where  suffix  o  denotes  initial  values  as  usual,  and  oscillations  of  amplitude 
less  than  10  km  have  been  ignored  (because  the  only  orbits  available  for  these 
satellites  generally  have  an  accuracy  in  perigee  height  not  better  than  10  km). 
If  we  write  the  last  three  terms  in  equation  (33)  as  Ay^j  ,  Ay^  ^p3  * 

their  variations  are  as  shown  m  Fig  26. 


However,-  y^  is  not  the  right  value  to  use  m  equation  (32):  we  need  the 
value  of  yp  that  correctly  gives  the  mean  density  at  perigee  over  a  cycle  of 
the  oscillation.  This  height,  y1  say,  is  always  someihat  less  than  yp 
wc  may  write 


and 


y 


1 


y  -  Ay  , 
P  P* 


(34) 


where  Ay^  is  as  given  in  Fig  27.. 
taking  the  perigee  huight  as  yl  . 


The  lifetime  can  now  be  read  from  Fig  25,. 


If  the  inclination  differs  from  25  ,  the  three  ntaocncal  values  m 

equation  (33)  change.  If  we  call  then  Aj  ,  and  A^  we  find  that  Aj 

decreases  from  43  kn  at  i  *  10  to  28  km  at  i  »  30  ;  in  contrast,-  A^  increases 

from  8  km  at  i  »  10  to  33  km  at  i  ■  ?0  ;  and  A^  increases  from  15  km  at 

l  »  10  to  26  k i  at  i  *  30  .  The  calculation  of  y  should  be  modified 

p 

accordingly  if  i  "*■  25  , 


The  method  is  valid  if  the  liletioe  extends  to  several  cycles  of  the 
perturbation,  that  is,  if  L  exceeds  about  3  years.  This  implies  y*  >  150  km 
2 

from  Tig  25,  if  m/S  «  100  kg/m  .  The  results  obtained  witn  the  method  should  be 
valid  to  within  a  factor  of  1.5,  subject  to  the  usual  pxoviso  about  solar 
activity. 
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10.3  Molniya-type  orbits  (iQ  *  65°  or  63°,  eQ  0.74) 

The  launches  of  Molniya  communications  satellites  with  12-hour  orbital 
periods  and  initial  inclinations  near  65°  began  in  1965,  and  31  such  satellites 
were  launched,  with  initial  eccentricities  near  0.74,  accompanied  by  31  separated 
rockets  m  similar  but  slightly  different  orbits.  Then  in  1973  there  was  a 
change  of  inclination  to  63°,  which  has  been  retained  for  all  the  82  launches  in 
the  years  1974-1986. 

The  perigee  height  of  a  Molniya  increases  from  about  500  km  initially  to 
more  than  1000  kia,  and  then  descends  some  years  later.  For  the  Molniyas  at 
65°  inclination,  an  approximate  analytical  solution  can  readily  be  devised  and 
the  results  are  recorded  here  m  case  launches  at  65°  inclination  are  restarted. 
For  the  63°  Molniyas,  the  lifetimes  are  touch  longer  and  less  amenable  to 
analytical  evaluation,  although  an  approximate  indication  can  be  given. 

For  the  65°  Molniyas,  numerical  integration  with  the  aid  of  the  PROD 
computer  program  shows  that  the  inclination  nearly  always  stays  between  64° 
and  66°,  and  i  may  be  taken  as  65°  with  an  error  in  sin  i  of  less  than  0.7 7,. 
The  changes  m  eccentricity  due  to  lunisclar  perturbations^6  can  then  be 
expressed  algebraically.  If  only  the  most  important  terms  are  included,  it  can 
be  shown6  that  Che  perigee  height  is  driven  down  to  100  km,  and  the  lifetime 
ends,  when 


0.566  +  20;'  -  0.31  cos  0 

_  o 

1  -  0.36  cos  0 


~  cos  2u  ±  0(0.04) 


(35) 


where  c'  is  ar  empirical  parameter  to  cater  *or  the  departure  of  the  initial 
eccentricity  from  0.74:  it  is  found  that  -0.005  <  c*  <  0.007.  If  we  exclude  the 
first  launch  as  non-standard,  29  of  the  remaining  30  M-'iniya  launches  at  65° 
inclination  had  *  285°  (with  one  at  284°).  Gracefully  accepting  this  devas¬ 
tating  demonstration  of  launching  accuracy,  we  may  take  •  285°  as  standard 
and  construct  a  diagram.  Fig  28,  showing  the  permitted  theoretical  lifetime 

domains  as  a  function  of  il  .  The  actual  decay  dates  of  29  of  the  30  satellites 
o 

are  also  shown:  all  are  either  inside  or  close  to  the  theoretical  domains. 
According  to  the  theory,  the  lifetime  of  a  65°  Molniya  wil*  not  exceed  7  years. 
The  Molniya  satellites  are  of  three  types  and  are  numbered  accordingly  m  Fig  28: 
thus  Molniya  1-12  is  the  12tn  launch  of  type  l.  In  Fig  28,  Molniya  1-23  has  been 
used  as  an  example  to  show  the  variation  of  lifetime  with  given  by  PROD. 
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The  theory  ignores  drag,  which  is  usually  important  only  in  the  last  month 
of  the  life.  But  the  drag  at  that  tine  can  sometimes  extend  the  life  by  a  few 
months,  because  the  drag  reduces  e  and  therefore  reduces  the  rate  at  which  the 
perigee  height  is  being  driven  down  by  lunisolar  perturbations. 

For  the  Molniyas  with  initial  inclinations  of  63°, ,  numerical  integration 
shows  that  the  inclination  may  subsequently  varv  between  62.8°  and  65.5°,  so 
that  w  varies  in  a  complex  manner,  usually  increasing  when  i  <  63.4°  and 
decreasing  when  i  >  63.4°.  It  is  possible  for  <n  to  become  trapped  near  2/0° 
for  several  years so  that  -cos  2u  remains  near  1.0  and  may  then  always  be 
greater  thar  the  left-hand  side  of  equation  OS),-  implying  an  indefinitely  long 
lifetime. 

The  results  for  63°  Molniyas  are  therefore  much  lesM  definite  and  nay  be 
summarized  as  follows**: 

(i)  Early  decay,  after  2  to  •*  years  in  orbit,  .s  possible  il 

135°  <  <  225°,:  but  it  is  unlikely  land  has  not  yet  happened,  except  for 

a  rocket) 

(li)  If  90°  <  <  240°,  there  is  a  good  uhance  (>502)  of  decay  between 

9  and  14  years  after  launch  . 

\ i i i )  If  -120  <0  <  90°,  Jeiav  is  1  ikely  between  1 2  ana  20  year s  after  launch, 

o  q 

unless  <ii  1  ibr«u  s  ibout  270  .  Il  libration  does  occur  (and  it  is  lik‘,|\ 
only  when  -60°  '<  no  upper  limit  can  he  specified  for  the 

lifetime. 

rig  29  shows  thes#  lue.irx.  domains  gr  tphic  illv.  v.ih  all  available  actual 
lifetimes  and  others  pr«*di<  led  b\  .inixi  u  i!  integration  with  I’ROI).  (The  .<*,  tual 
lifetimes  ire  ill  less  tine  15  \carx,  btiluse  tV  *>1°  '<oIiii\as  did  not  start 
until  1971.)  In  view  of  the  gross  approximit  ut.s  ,  .t  is  surprisin,'  that 
all  the  points  are  v* thin  Hu  dunams  indicated  m  Fig  23:  indeed  it  acens  that 
the  domains  might  be  made  smaller. 

foe  re  are  other  ,  <eei.il  high-e<  centru  it'  orbits  hut  they  ir<  (a)  m  t  -.o 
heavily  populated  is  the  transfer  and  Mol.iiya-typt  orbits,  aid  (b)  mostly  not  so 
amenable  to  treatment  by  approximate  methods,  because  the  inclination  often 
varies  rather  widely. 


Note  that  the  upper  limit  for  0o  has  been  revised,  from  270°  m  Ref  6  to 
240  here,  in  the  light  ot  experience. 
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11  CONCLUSIONS 

Anyone  who  has  struggled  through  to  the  end  of  this  Report  will  emerge  with 
some  idea  of  the  difficulties  and  complexities  of  predicting  the  decay  date  of  a 
satellite.  Accurate  lifetime  prediction  is  one  of  the  most  incractable  problems 
of  orbital  mechanics. 

It  cannot  be  said  too  often  that  there  is  no  glory  to  be  gained  from  life¬ 
time  prediction.  If  you  do  succeed  in  making  an  accurate  prediction,  most  people 
think  it  must  be  easy  and  the  rest  will  say  it  is  just  luck.  If  your  prediction 
is  wrong,  you  can  expect  sneers  from  the  ignorant  and  the  semi-learned. 

We  favour  the  grpahical  methods  presented  here,  because  they  are  simple 
and  of  adequate  accuracy  and  comprehensiveness.  Others  may  prefer  more  sophisti¬ 
cated  methods,  such  as  numerical  integration;  but  their  predictions  are  unlikely 
to  b,_  more  accurate,  because  the  errors  arise  chiefly  from  the  future  variations 
in  air  density,  which  are  uncertain. 


i. 
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Appendix 

ORBITAL  LIFETIME  IN  AN  ATMOSPHERE  WITH  DAY-TO-NIGHT  DENSITY  VARIATION 

We  assume,  as  in  Ref  22,  that  the  density  p  at  distance  r  from  the 
Earth's  centre  is  given  in  terms  of  the  bulge-perigee  angle  $  by 


p  *  p'(l  +  r’  cos  $)  expj-  (r  -  r')/H|  , 


(A.  1 ) 


t 


where  F  is  a  constant,  equal  to  (pu  -  p  >/(pM  +  p  )  in  the  notation  of 

n  m  H  m 

section  7;  $  is  as  defined  in  Ref  22;  and  the  constant  o'  is  the  density  at 

distance  r'  from  the  Earth's  centre  for  $  *  90°,  where  the  value  of  r'  is 
open  to  choice  (it  would  normally  be  taken  as  the  initial  perigee  distance  r^Q 
but  this  is  not  compulsory) . 

If  T  (»l/n)  is  the  orbital  period,  we  have  T  »  3Aa/(2a),  where  Aa  is  the 
change  in  a  during  one  revolution,  and  hence,  from  equation  (35)  of  Ref  22, 


-  ter 


a6p  {l  +  2e  +  - 


2  0.1  F  F*  \l 
V  ’  2  •  2>  2z)l 


(A. 2) 


for  Phase  1,  that  is,  fot  3H/a  <  e  <  0.2.  Here  6  is  the  area/mass  parameter 
defined  m  section  9.  In  view  of  the  other  uncertainties  m  the  analysis,  we 
shall  neglect  terms  le3S  than  0.1  in  equation  (A.2).  Since  the  highest  value  of 
F  attained11  is  0.7,'  we  lose  little  by  assuming  that  F  <  0.6,  so  that 
F  cos  $p/(2z)  <0.1.  Also  l/(8z)  <  0.05.  Thus  equation  (A.2)  reduces  to 

i  -  -  3(1/21)*  »5pp{l  *  2e  *  0(e2,  0.1/z,  iF/z)}  .  (A. 3) 

Following  the  method  of  Ref  4,  we  write  down  the  corresponding  equation  for 
the  unperturbed  orbit  (with  F  *  0),  denoting  values  on  that  orbit  by  suffix  u  : 

K  -  -  3<*/2*u>!  w{'  *  2% +  0('2’  •  (A-4) 

We  may  drop  the  suffix  u  on  z  ,  a  and  e  ,  because  the  maximum  change  in 
(r^Q  -  rp)/H  due  to  the  day-to-night  variation  is,  from  equation  (72)  of  Ref  22, 

of  order  ^-F  cos  $  ,  which  is  less  than  0.1  for  F  <  0.6;  and  it  then  follows 

.  p 

from  equation  (15)  of  Ref  4  that  zy  ,  and  e^  are  not  significantly  differ¬ 
ent  from  z  ,  a  and  e  .  Also,  from  equation  (A.l)  with  F  -  0, 

t  The  symbols  r^  and  pq  in  Ref  22  have  been  replaced  by  r’  and  p'  ,  to 
avoid  confusion  over  the  suffix  o  . 
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.  Hv  -  y/H} 


(A. 5) 


Thus  equation  (A. 4)  becomes 


-He) 


3Sppou{‘  +  2c  *■  °^2 


°-*«}  exp{(rpo  '  rp,/H} 


.....  (A. 6) 

•  4 

We  require  that  T  «  initially  :  comparing  (A. 3)  and  (A, 6),  we  see  that  this 

requires  p  to  be  defined  as  equal  to  P  ^  •  As  Ref  4  shows,  F(dn)  is 
pou  .  .  po 

given  by  the  mean  value  of  T^/T  during  the  remaining  life  of  the  satellite. 
Thus,  dividing  (A. 6)  by  (A. 3)  and  noting  that  p  is  constant,  we  have 


F(dn) 


/■A. 

ppo 

\T  f  meanjp 

P  H(rP  -  rpo)/H}] 

(A.  7) 


The  denominator  here  is  the  mean  value  of  the  perigee  density  at  the  initial 
perigee  height  (but  of  course  at  the  local  value  of  ^  ) »  and  represents  the 
quantity  denoted  by  in  equations  (8)  and  (10).  Thus  F(dn)  is  as  given 

in  equation  (15).  Since  equation  (15)  is  the  same  as  equation  (8),  and  the 
solar-cycle  variation  of  density,  given  by  equation  (9),  has  the  same  form  as 
the  da> -to-night  variation,  given  by  equation  (13),  the  results  obtained  for 
F(sc)  apply  also  for  F(dn)  ,  if  0.03  <  e  <  0.2. 

The  only  formal  limitation  that  has  accrued  during  'he  analysis  is  that  c 

should  be  'Phase  1*,  approximately  G.03  <  e  <  0.2.  However,  there  is  another 

limitation:  the  mode]  for  density,  (A. I),  will  not  always  be  adequate  if  the 

satellite  decays  before  completes  one  revolution,  because  there  may  be  a 

great  decrease  in  r  ,  so  that  the  values  of  F  and  H  ought  to  be  con^ider- 
F 

ably  reduced.  This  caution  is  not  needed  if  F  <  0.2,  as  happens  if  the  perigee 
height  is  less  than  200  km. 

The  result  (A. 7)  obtained  for  Phase  *  applies  a  fortiori  for  e  >  C.2 
because  the  terms  in  \iz  become  quite  negligible,  and  the  terms  in  c  and 
m  equation  (A. 2)  can  be  grouped^  as  (1  ♦  e)^^(l  -  e)  ^  without  requiting  the 
assumption  that  e  is  small. 

For  e  <  0.03,  the  analysis  takes  a  slightly  different  form.  From 
equation  (31)  of  Ref  22,  we  have 
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T  -  -  3*a5p 1  ^Iq  ♦  Ij  F  cos  +  0(e,  eF;j  expj(r'  -  a)/tlj 


(A.  8! 


From  equation  (A.l), 

Opo  ■  p,(l  +  F  CM  V  expi<r’  '  rpo)/H}  •  <A-9 

Dividing  equation  (A. 8)  by  equation  (A. 9)  gives 


T  *  -  3xa<5p 


In  ♦  F  cos  *  *  0(e,  Fe) /  ir  -  r 


1  ♦  F  cos  + 


X2 _ 2.  _ 

H 


(A.  JO 


On  the  unperturbed  orbit,  equation  (A. 10)  takes  the  form 


*u  ’  -  ^“pou^o  +  °<*>} 


r  -  r 
_ L-  ; 


(A.  11 


where  we  have  again  assumed  that  z  and  a  do  not  differ  significantly  from 
their  unperturbed  values,  relying  for  this  assumption  on  equation  (130)  of 
Ref  22,  which  is  similar  in  its  main  term  to  equation  (72)  of  Ref  22.  (There 
may  be  problems  with  this  assumption  as  z  “*  0  :  but  then  the  effect  of  the  day- 
to-night  variation  tends  to  zero.) 

Dividing  equation  (A. II)  by  equation  (A. 10),  we  have 


T  /P  \  l  +  F  cos  4  .  y 

f  '  +  0<e’ Fe>}  •  (*’12! 


The  value  of  Ppou  is  defined  by  the  requirement  that  T  ■  Ty  initially.  Thus 


t 

u 

T 


1  *  <W  F  ”  Sc 

I  +  (X./In)  F  cos  *  {'  0(e’  Fe)f  ’ 


(A.  13) 


where  I  is  the  initial  value  of  I  . 
no  n 

To  obtain  F(dn)  for  e  <  0.03,  we  integrate  T  /T  over  the  remaining 
lifetime  L  ,  as  in  Ref  4,  so  that 
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F(dn)  •  ~  /  [l  '  (1,/y  F  V  *  (II0/IO0)  F  cos  vl  dt  • 

0 

(A. 14) 

i£  we  assume  that  F  ia  small  en.  igh  for  terms  in  E2  to  be  ignored.  In 
attempting  to  integrate  equation  (A.  14) ,  the  obvious  first  approximate  .s  to 
assume  that  !,/!<>  rematns  constant  at  tts  initial  value,  (the  decrease  of 
I  /I  with  time ,  shown  in  Fig  30  for  so  -  3,  assuming  z!zq  -  (1  -  t/»  .  is 

Ur/slow  until  near  the  end  of  the  life.)  Then  assuming  that  ?p  vanes 
linearly  vith  time,  with  *p  •  2a/?  ,  equation  (A.14)  gives 

F(dn)  -  I  *  f[ros  *  -  ~l  {Sln(>  *  ¥}'  0] 

00  L 

(A. 15) 


The  corresponding  equation  for  0.03  <  e  .<  0.2  is  obtained  on  dividing 
equation  (A.6)  by  equation  (A.  3),  to  give 

Tu/T  *  <V'V  HV  '  r?)/H|  ’  1 

since  P  -  0  Phase  1,  as  already  noted  after  equation  (A.6).  From 

>P0U  po 
equation  (A.l)* 


J2£ 


■  *  P  tos  »po 

1  +  F  cos  ♦p 


exp{(rp 


o>/H} 


(A. 17) 


From  equations  (A. 16)  and  (A. 17), 


T 

u 

T 


1  ♦  F  cos  »pQ 
1  +  F  cos 


(A. 1 8) 


comparison  with  equation  (A. 13)  shows  that  the  result  for  0.03  <  e  <  0.2  is  the 

same  as  for  e  <  0.03,  except  that  I,/I&  is  replaced  by  1 .  Thus  if  F  terms 

are  ignored,  F(dn)  for  0.03  <  e  <  0.2  is  given  by  equation  (A. 15)  with 


030 


030 


Consequently,  when  <  3  ,  the  'strength'  of  F(dn),  ie  its  departure 
from  1.0,  needs  to  be  scaled  down  by  a  factor  Ij(z  )/Iq (z  )  and  is  given  by 
equation  (17).  Values  of  Ij/IQ  appear  in  Table  A.l  at  the  end  of  this  Appendix. 

In  integrating  equation  (A.  14)  for  z^  <  3,  we  have  assumed  that  Ij/Iq 
constant.  This  is  an  acceptable  approximation  if  zq  is  near  3,  but  over¬ 
estimates  the  'strength*  of  F(dn)  ,  particularly  for  z^  <  1.  An  opposite 
extreme  is  to  assume  that  Ij/Iq  decreases  linearly  to  zero,  and  Fig  30  shows 
that  the  truth  lies  between  these  approximations.  Taking 


(,  .i) 

Inn  V 


in  equation  (A.  14),  we  find 


— —  F  I" cos  $  +  sin  $  +  -SL.  Losfa  ♦  -  cos  * 

T00  l  po  2rL  po  4.2L2  '  '  p°  f  '  p0’J 


Two  of  the  terms  within  square  brackets  in  equation  (A. 15)  are  the  same  as  in 
equation  (A. 19),  and  the  third  differs  from  the  remaining  terms  in  equation  (A. 19) 
by  less  than  0.1  if  L/P  >  J  and  F  <  0.3.  Thus  the  two  approximations  give 
effectively  the  same  result. 


In  deriving  the  results  for  e  <  0.03,  we  have  ignored  the  term 
{(F  Iq/Ij)  cos  $  {2  ,  which  integrates  to  give 


4©  <'*•<&»  • 


This  term  is  less  than  0.1  if  F  <  0.55  for  z  ■  3.  or  if  P  <  0.64  for  z  *  2,  so 

o’  o’ 

its  neglect  is  in  close  accord  with  our  previous  assumption  thar  F  <  0.6. 
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Accurate  prediction  of  the  lifetimes  of  satellites  decaying  naturally  under  the 
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For  use  in  the  RAE  Table  of  Satellites  we  have  since  1958  made  lifetime  esti¬ 
mates  for  all  satellites,  rockets  and  other  major  components.  More  than  10000  such 
ipredictions  have  been  made,  nearly  all  being  calculated  from  the  current  observed 
(decay  rate.  This  Report  summarizes  the  methods  used  in  the  past  and  gives  the 
results  of  new  orbital  theory  developed  to  take  better  account  of  the  effects  of 
atmospheric  ob lateness  and  odd  zonal  harmonics  in  the  geopotential.  The  Report  also 
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to-night  variations  in  air  density,  and  includes  discussion  of  other  sources  of 
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